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Abstract FWHM calibration is one of the essential procedures in gamma-ray spectra analysis. Traditional method is 


time consumption and poor flexibility because of extra measurement on calibration source. This paper introduces an 


implicit calibration method based on symmetric zero-area conversion. The single peaks can be differentiated from 


doublets. Meanwhile, the extracted FWHM values can be applied in the fitting process. Various forms of calibration 


equations are introduced and compared with practical spectrum data under this method. 
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1 Introduction 


The particular measure of the peak widths in gamma- 
FWHM (Full width at half 
maximum), indicates the resolution of the gamma 


ray spectrometry, 


spectroscopy system. The FWHM expressed in keV 
has a relationship with energy of incident gamma-ray. 

It is important to know the FWHM of peaks as 
an energy function when analyzing gamma spectra, 
especially in peak searching and fitting!" Thus, an 
accurate FWHM calibration needs be performed as the 
basis of spectrum processing. 


2 Physical essence of FWHM 


Gamma-ray spectrum is actually the energy response 
of the detection system to the gamma-rays of 
particular energies. Peaks in the spectrum are spread 
over several channels that can be characterized as a 
Gaussian distribution, whose standard deviation can be 
predicted by”! 
O=0 1+0 pto cto p (1) 
o is the overall uncertainty of the peak in the 
spectrum in energy form. 
aris the intrinsic width of nuclear energy levels, 
which is extremely smaller than o, and can be 


ignored"), 
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op is the statistical fluctuation of the production 
of electron-hole pairs (n) in the detector. The number 
of electron-hole pairs created by a gamma-ray of 
energy (E) subjects to Fano distribution (PP! On 
defining "W" as the average energy consumption to 
create a pair, the expected uncertainty of n is 
on=(Fxn)” = (FXE m]. 

oc is the uncertainty in collecting the charge by 
the detector caused by trapping effect. Due to the 
complexity effect, there is no a simple way to express 
the uncertainty as an energy function®?. However, a 
linear relationship of oc=cE appears to be satisfactory, 
but not theoretical reason"). 

of, the contribution from electronic noise in 
pulse processing system, is performed as a constant to 
the gamma-rays with any energy. 


3 FWHM calibration equations”! 


There are mainly three independent effects in 

influencing the resolution of gamma-ray spectrum. 

Using the energy form AF to express the FWHM of a 

Gaussian-shaped peak in the spectrum, we can give 
AEX2.3550" 


Or AE =(P E) cE) He? (2) 


where p, c, and e are constant related to fluctuation, 
collection, and electronic noise, respectively. 
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The square root for both side of Eq.(2) is taken, 
And 


AE=\pE+CE’ +e = a +a? E+E (3) 
Debertin and Helmer"! introduced Eq.(4) by 
ignoring the factor of incomplete charge collection. 


AE =; pE +e = a? +a E (4) 


In addition, Eq.(5) with simpler expression can 


Table 1 Results of fitting FWHM data to different functions!” 


be found in the commercial programs of Genie 2000", 
AE =a, + aNE (5) 


Gilmore conducted a series of experiments to 
compare the fitting FWHM calibrations with different 
formulae”! (Table 1). The best fit is from the square 
root quadratic function, and a simple quadratic fit. The 
Genie 2000 fitting is worse than the linear equation. 


Fitting type Formulae RMS differences 
Linear agta,E 0.020 

Quadratic agta, E+ta E? 0.0065 

Genie 2000 ayta,E°° 0.055 

Debertin and Helmer (ao +a rE 0.026 

Square root quadratic (aota, E+aE°) 0.0057 


4 Traditional FWHM calibration 


Similar to energy calibration, traditional FWHM 
calibration method is accomplished by measuring 
gamma-rays of precisely known energy and estimating 
the relationship between the peak width and energy". 
The three ways can provide proper gamma-ray 
emission source. (1) Several kinds of nuclides with 
single or discrete energies are chosen and mixed 
together; (2) One specific nuclide emitting a variety of 
energies is chosen, covering the entire range of the 


2h. The known nuclides in the 


spectrum, such as 
test sample can be used to perform the calibration. 
Explicit method can be used to fine calibration; 
but its calibration source needs be additionally 
measured by time consumption. Especially, when 
changing experiment condition, the spectroscopy 
should be 


measurement except applying for nuclides in any 


recalibrated by employing a new 
situation. In this paper, an implicit method for FWHM 
calibration is developed by self-information gained 
from the target spectrum. 


5 Implicit FWHM calibration method 


5.1 Method introduction 


The Ge(Li) detector has high energy resolution, and 
makes a great number of gamma peaks with enough 
distance each other, separating single peaks easily. 

two calibration 


Implicit method contains 


procedures. First, single peaks are examined and 


extracted from the spectrum. Second, the width of 
each singlet is calculated. Together with the peak 
position, the calibration curve is obtained by fitting 
these data. 

Comparing all kinds 


implicit 


of peak 
FWHM 
symmetric zero-area conversion method (SZA) has 


searching 
algorithms with calibration, 


advantages in low noise sensitivity, fine multi-peaks 


resolving and weak peak identifying abilities!) 


5.2 Method description 


The spectrum data are convoluted by the SZA method 
using a zero-area window. The linear base tends to 
zero after convolution while maximum value is at the 
positions of gamma peaks. The conversion process is 
described as Eq.(6). 


= 6) 


where y; and y; are the spectrum data before and after 
conversion, and C; is the window coefficient. 

Among various forms of window functions, 
Gaussian second derivative has the best property’, its 
coefficients can be expressed as Eq.(7). 


WA)=C,=—? -ode -70o (7) 
om 2 


where, o,, is the standard deviation from the original 
Gaussian distribution G,,(/). 
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Figure 1 shows a Gaussian second derivative 


formed SZA window function with o,,=4. 
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Fig.1 Gaussian second derivative formed SZA window 


function. 


And, gamma peaks in the spectrum can be 
treated as Gaussian shape approximately. 


20, 


G, (i) =y, = meon — py | (8) 


where, H, p and o, are their amplitude, position and 
standard deviation, respectively. 


Due to the convolution naturel”, a Gaussian 
peak conversion operated with the window function 


second 


of Gaussian 


maintains the appearance 


derivative, as shown in Fig.2. 


Fig.2 Convolution between Gaussian peak and window 
function. 


RO=G (WG) =C) G7)" 
=[G,0)* G7) = Geonv(i)" (9) 
Geony(i) is the convolution of gamma peak and 
the original Gaussian distribution of the window 
which has a Gaussian form of d¢ony= oe) 
The distance between two zero-crossing points 
beside the conversion maximum can be obtained from 
the Gaussian second derivation. 


2 2 \0.5 
D=26conv=2(6 wto z) 


(10) 


The o, of the gamma peak is deduced from 
intercept D and cw, subsequently, the peak FWHM is 
expressed by Eq.(11). 


FWHM = 2.3550, =2.355 (D° / 4)- o} 


The Eq.(11) is only valid for single peaks, and 


(11) 


the calculated FWHM is definitely a constant for every 
pair of D and o,. While for multi-peaks, the intercept 
D relates to the intervals and ratios between the 
component peaks, in other words, the FWHM obtained 
from Eq.(11) is meaningless. 

Consequently, the same spectrum can be 
converted using several windows with different ow. 
When the FWHM value tends to be invariable, the 
corresponding peak is considered as a singlet, thus 
distinguishing a set of multi-peaks due to their invalid 
FWHM. Also, the FWHM of each single peak can be 
gained, implementing the FWHM calibration process. 


5.3 Method modification 


Applying to practical spectra, a group of simulation 
data set up is used to test the calibration method. 

The o,, of 8 window functions equal to 1 to 8 is 
convoluted by o, of 4 individual Gaussian peaks equal 
to 1, 2, 4, and 8. The calculated and theoretical FWHM 
values are compared (Table 2, Fig.3). 
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Fig.3 Error between the calculated and theoretical FWHM 


values under different windows. 


Table 2 Difference between calculated FWHM and theoretical value 


Error between calculated and theoretical FWHM values 


Theoretical 
K FWHM ov=l Oy=2 Oy=3 O„=4 O, =5 O„=6 Oy=7 o,=8 Mean Std 
1 2.355 0.321 0.279 0.204 0.152 0.116 0.0881 0.0650 0.0445 0.159 0.101 
2 4.710 0.146 0.0542 0.124 0.142 0.139 0.127 0.114 0.0992 0.118 0.0301 
4 9.419 0.0405 0.0700 —0.00681 0.0527 0.0646 0.0442 0.00736 -1.71e4 0.0340 0.0299 
8 18.84 0.00967 0.0258 0.0263 —0.00798 0.0177 —0.0202 0.0114 0.0102 0.00912 0.0161 
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To reduce the difference in Eq(1l), a 


modification factor is introduced to Eq.(12). 


FWHM =2355 (D° /4)— 0} -(a/ Jo,) (12) 


The constant (a) is determined by applying an 
iteration method to optimize the average error under 
every window function for each Gaussian peak. 

Figure 4 shows the result using Eq.(12). Under 
every window with different o,,, the mean of residual 
errors approaches to zero, and the fluctuation is 
suppressed at some levels. 

In practical application, the o,, values should be 
selected. The SZA method can separate multipeaks 
with small o,,, but is sensitive to the noise. On the 
other hand, the reliable window of large o, is 
inadequate in multipeak processing. For ow, a value 
approximating to o, suggests obtaining the optimal 
results. Therefore, o, can be determined within a 
rough range related to FWHM varying in the spectrum. 
Also, this method is compared with calculation 
process. The average background level above and 
below the peak is subtracted from the estimated peak 
height to get the expected count. The two positions of 
half height on the peak are obtained by linear 
interpolation, and the distance is exactly the FWHM. 


Table 3 Compare of direct process and method in this paper 


The resultant reliability depends on the uncertainty of 
peak height and background level. 

Quadratic formed background and noise with 
constant SNR are added to a Gaussian peak by o,=4 to 
simulate an actual gamma peak. The 100 sets of 
simulation data with random background and noise are 
processed by the two methods (Table 3). 

The FWHM Eq.(12) 
approximates to the theoretical value. For the direct 


calculated by 


approach, the result is relatively fluctuant due to the 
noise. Applied first to reduce the standard deviation, a 
smoothing process impacts the FWHM value to a 


certain extent. 
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Fig.4 Error between calculated and theoretical FWHM values 
Eq.(12). 


Error of 100 experiments 


Method Theoretical FWHM 

Mean Std 
Direct-without smoothing —0.0296 0.252 
Direct-with smoothing 9.419 —0.0405 0.221 
This paper 0.0167 0.170 


6 Method verification 


A gamma-ray spectrum from sand sample is measured 
by HPGe detector, and used to test the implicit FWHM 
calibration method (Fig.5). 
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Fig.5 Practical gamma-ray spectrum. 


The 3 windows with o,, equal to 2, 3, and 4 are 
employed to convert the spectrum and calculate the 
FWHM corresponding to each peak, the three sets of 
FWHM deviation is applied to pick out single peaks, 
and their mean value is used to the fitting procedure. 

HYPERMET, a mature developed y spectrum 
analysis code, and Genie 2000 program are employed 
to validate the in Table 4. The 
HYPERMET using a region fitting method to get 
accurate peak position and FWHM can be used as the 


calculation 


standard value, while Genie 2000 applies a direct 
process to calculate FWHM. Together with FWHM 
values, the relative errors comparing to HYPERMET 
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are calculated for the other two processes, as shown in both HYPERMET and Genie program need a rough 
Table 5. The relative errors within +3% are obviously FWHM calibration as input to obtain proper results, 
better than Genie 2000 except one point. In addition, while the novel method is fully automatic. 


Table 4 Calculated FWHM and position of each single peak 


a Calculated FWHM 
Peak position 
Oy=2 Ow=3 Oy=4 Mean Std 
323 2.729 2.917 2.833 2.826 0.0946 
470 2.770 2.835 2.773 2.793 0.0368 
631 2.990 3.111 3.110 3.070 0.0699 
812 3.236 3.254 3.267 3.253 0.0153 
1116 3.410 3.414 3.398 3.407 0.00824 
1407 3.635 3.653 3.592 3.627 0.0312 
1531 3.782 3.722 3.748 3.751 0.0301 
1918 3.976 4.025 4.065 4.022 0.0443 
2197 4.227 4.217 4.220 4.222 0.00494 
2861 4.600 4.616 4.634 4.616 0.0169 
4127 5.446 5.538 5.501 5.495 0.0460 
4335 5.580 5.642 5.583 5.602 0.0347 
4595 5.754 5.847 5.703 5.768 0.0728 
6851 7.049 7.131 7.229 7.137 0.0901 
Table 5 The FWHM calculated by different methods 
fe Calculated FWHM 
Peak position i : ; 
HYPERMET Genie Relative error /% This paper Relative error / % 
323 2.751 3.024 9.939 2.826 2.726 
470 2.905 2.945 1.371 2.793 -3.855 
631 3.131 3.210 2.527 3.070 -1.948 
812 3.190 3:237 1.463 3.253 1.975 
1116 3.368 3.396 0.8266 3.407 1.158 
1407 3.660 3.688 0.7560 3.627 -0.9016 
1531 3.685 3.714 0.7924 3.751 1.791 
1918 4.025 4.245 5.461 4.022 -0.07453 
2197 4.190 4.218 0.6747 4.222 0.7637 
2861 4.588 4.457 -2.854 4.616 0.6103 
4127 5.446 5.651 3.762 5.495 0.8997 
4335 5.550 5.571 0.3874 5.602 0.9406 
4595 5.727 6.261 9.326 5.768 0.7159 
6851 7.065 7.137 1.007 7.137 1.013 
Table 6 and Fig.6 show the compared position- 5 
FWHM data from Table 4 fitted to the 5 equations in 
Section 3. The square root quadratic function has the 77 ] 
best effect, and the linear fit performs better than 6l | 
Genie 2000 function, this is consistent with the = 
conclusion in Section 3. z > 
Table 6 Compare of different equations fitting FWHM data 47 ea ey | 
“Calibration equation RMS differences 3} ~~ Gene 2000 | 
Square root quadratic 0.0235 f R oak 
Quadratic 0.0285 * 1000 2000 3000 4000 5000 6000 7000 
Linear 0.0774 Channel 
Debertin and Helmer 0.162 
Genie 2000 0.433 Fig.6 Fitting curves of the 5 calibration equations. 
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The performance of above 3 typical equations calculation process first using simulation data. 
is shown in Fig.7. Ultimately, this method is applied to a practical 
0.4 r r z gamma-ray spectrum, and the calibration result is 


—@-— Linear 


E E ae E pons proved to be good. 
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